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Abstract 

We calculate the 0{as) SUSY-QCD corrections to the widths of stop and sbottom de- 
cays into Higgs bosons within the Minimal Supersymmetric Standard Model. We give the 
complete analytical formulae paying particular attention to the on-shell renormalization 
of the soft SUSY-breaking parameters. We also perform a detailed numerical analysis of 
both stop and sbottom decays into all Higgs bosons h^, H^, and H^. We find that 
the SUSY-QCD corrections are significant, mostly negative and of the order of a few ten 
percent. 



1 Introduction 

Supersymmetry (SUSY) requires the existence of two scalar partners g^, and qj^ to each quark 
q. In the case of the scalar partners of the top quark one expects a large mixing between ti 
and tfi due to the large top quark mass The mixing of and bji may also be substantial 
if tan/? = V2/V1 is large (where Vi and V2 are the vacuum expectation values of the two Higgs 
doublets). Strong mixing induces large mass differences between the lighter mass eigenstate qi 
and the heavier one q2, q = t or b. This implies in general a very complex decay pattern of 
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the heavier states. In addition to the "conventional" decays into neutralinos, charginos, and 
gluinos {ij = 1,2; /c = 1, . . . ,4) 

ii ^txl, bXj, hi bxl, txj, (1) 

ii^tg, bi-^bg, (2) 
decays into vector bosons and Higgs particles can become kinematically possible = 1,2): 

i2^hZ^ i,^b,W+, h^hZ^, h^i.W-, (3) 
t2 -^ti(/i°,if°,v4°), ii-^bjH+, b2^k{h\H\A^), bi~^i,H-. (4) 

All these squark decays were first discussed at tree-level in within the Minimal Supersym- 
metric Standard Model (MSSM) [Q. A recent, more complete and systematic analysis of these 
decays at tree-level in [^] revealed that the bosonic decays of Eqs. and can be dominant 
in a wide range of the MSSM parameters due to the large Yukawa couplings and mixings of t 
and b. This could have an important impact on the search for t2 and 62 and the determination of 
the MSSM parameters at future colliders. Therefore it is important to study how SUSY-QCD 
corrections affect this tree-level result. 

Within the last years SUSY-QCD corrections to a variety of processes were calculated. For the 
decays of Eq. (|l]) this was done in P, ||, 01 > for the decays Eq. (||) in § , and for the decays 
Eq. (j^) in 0]. The SUSY-QCD corrections for the decays into Higgs bosons, Eq. (^, were briefly 
discussed in [|10[. The QCD corrections to the related Higgs boson decays {h^,H^,A^) q^q^ 



and q^q, were calculated in |]lT], |10[. A thorough study of the corrections to the decays 



Eq. (p, including a detailed numerical analysis, is, however, still missing in the literature. 

In this article we discuss the 0{as) SUSY-QCD corrections to the decay widths of Eq. (^) 
in the on-shell renormalization scheme within the MSSM. We give the complete formulas for 
these corrections and point out some subtleties which occur in the on-shell renormalization 
scheme. Whereas a numerical analysis was made only for ^2 — > + {h^, in '^^) we perform 



a detailed analysis on stop and sbottom decays into all Higgs bosons ^4°, and H^. 



2 Tree— level formulae and notation 

We first summarize the tree-level results and our notation. The squark mass matrix in the 
basis (g^, g^) is given by H 



with 



"^L = M? +m|cos2/5(J|i-egSin2^H/)+m2, (6) 
"^L = ^{u,D} + "^1 2/? Cq sin^ Ow + ml, (7) 
ttq = Aq — fi {cot P , tan P} (8) 
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for {up, down} type squarks. Mq jjj-, and At^b are soft SUSY-breaking parameters and is the 
Higgs mixing parameter in the superpotential. /| and are the third component of the weak 
isospin and the electric charge of the quark q. The squark mixing matrix TZ^ is 

/ cos 9 a 



sin dq 
sin 6q cos Oq 



(9) 



with < 0g < TT by convention. The weak eigenstates Qj^ and are related to their mass 
eigenstates and ^2 (with rriq^ < rriq^) by 



7^« 



(10) 



In the (g^, basis the squark interaction with Higgs bosons Hj^ = {h^, ^4°, H"^} can be 
written as {i, j = 1,2; k = 1 . . . 4; a and /3 flavor indices) 



r ryoi tt] ~P] ~a 

l-qqH — ^ijk'^k'ij Hi ■ 



The couphngs Gfjj^ are 



fig 

^ijk 



fit 



b 



(fc = 1,2,3), 



:ir 



:i2) 



(13) 



^, n. — 1 ... 4, being the respective couplings in the (g^, q^) basis as given in the Appendix. 



with Gf^, k 
Note that G^j^ 
is thus given by 



Gjj4 and Gj^g 



The tree-level width of the decay gf 



r°(gT 



where m,- 



[{x 
4 we have 



QjHk, Fig. la. 



(14) 



m-/3, rriH = rrij^ , and k(x, y, z) = [{x — y — — Ayz]^/"^. For k = 1, 2, 3 



we have of course a = [3 and i = 2, j = 1. For A; = 4 we have {qf,qj) = {ti,bj) or {bi,tj). In 
the following, we will omit flavor indices when possible (flavor = a if not given otherwise). 



3 SUSY-QCD corrections 

The 0{as) corrected decay amplitude in the on-shell renormalization scheme can be expressed 
as 

^ij^" = ^ijk + f^G-j]. + + (15) 

where the superscript v denotes the vertex correction (Fig. lb), w the squark wave-function 
correction (Fig. Ic), and c the counterterm due to the shift from the bare to the on-shell 
couplings. The 0{as) corrected decay width F is then given by 

F = r° + (5F(^) + (JF^"") + (5F(^) + STreal (16) 
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with 



ST^"^ = ^ ; ^ Re {G*., 5G ; } (a = t;, c) (17) 

and 5r^eai the correction due to real gluon emission (Fig. le) which has to be included in order 



to cancel the infrared divergences. We use dimensional reduction [|T2[ as regularization scheme. 
Analogous calculations were performed for the crossed channels of Higgs decays into squarks 
in 
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3.1 Vertex corrections 



The vertex correction due to the gluon-squark-squark loop in Fig. lb is 



5G. 



= ^ G^jk [B,{ml 0, ml) + B,{m], 0, m]) - B^{ml^,ml m]) + 2 X Cq 



(18) 



with X = mf + m| — . Bq and Cq are the standard two- and three-point functions ||T3| 
for which we follow the conventions of [|T^. In this case, Cq = Co(m^, m\j^,m?-] A^, m^, m|). As 
usual, we introduce a gluon mass A to regularize the infrared divergence. 

The graph with the gluino-quark-quark loop in Fig. lb leads to 



5G. 



2U 



2 tts „ . c 

— - — rrig cosZtJq 

3 TT 



Bo{ml,ml ml) + 5o(m?^, m?, mj) + {4ml- J Co (19) 



for the decays into and {£ = 1,2), 
= -^^.^{m,sin2e, 

for the decay into A^, and 
2 a, 
3^ 



5o(m|,m?,m2) 



Boiml,mlml) + imi-ml)Co 
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Eo(m|,m?,mJ) + Bo{ml^,ml,ml) - m^Co] } , 



(20) 



2as ( [ , . 

- — i {mg<.Au + 



mqi3A22) s2 + {mqcA22 + mqiiAn) si 



T3 I 2 2 2\ 

Bo{mH+,mt,mf^) 



+ 
+ 
+ 



(m,a All - mgA2i) s2 + {mgc.A22 - m^Aia) 4 5o("^i , m^, m; 



"^5^21) s? + (m^^Aii - mqAu) s 



{mgfiA22 

m n {mlo. - ml + m?) {A22 + An si) 



Boim^.m^, m 



with Ar 



+ mg (77^+ - m^c - m^p) {A21 s" + A12 si 
- 2mg mqc mqp {A12 + A21 si) Cq | 

'^?rJ^jm the decay into a charged Higgs boson. 



(21) 



In Eqs. ([T9|) to ( pT|) Co = Co(mf , ""^1/^., m?, m^a, m^^)); are the Higgs couplings to quarks: 
CqqH = si h^qq + si H^qq + si A^gT^ q + H+ t{s\ Pr + 4 Pl) b + R-b (4 Pr + 4 P^) t (22) 
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with 



s = - 



^^^^ cos a, 4 



2 mw sin /3 



1 2 myi? cos /3 



sm a, 



So = - 



, ^ sin a, 4 = i^cot/5, 
, ^ COS a, s^. = i^tan/3 

2 miY cos p ' J 2 mi4/- ' 



^^COt/?, 



V2 



tan /5. 



(23) 



The vertex correction due to the four-squark interaction in Fig. lb is 

^{v,q) _ '^s ca C/3 



SG. 



ijk 



ra,m=l,2 



with 



cos — sin 26q 



sin 26*5 



cos 29 q 



(24) 



(25) 



3.2 Wave— function correction 



The wave-function correction is 



_ 1 



5Zt, + 5Z^^ 



where the are the squark wave-function renormahzation constants for They stem from 
the gluon, gluino, and squark loops of Fig. Ic Q and are given by: 



5^ir^ = -Re{E(f/)(m|j} , 5Z(f 



■7(9'^) 



Re{sgf)(m^-J} 

2 2 
-rrig 



n ^ n 



with S„„(m^) = 9S„„(p^)/9p^|p2=m2. The contribution due to gluon exchange is 

and those due to gluino exchange are 

2«, 

3 TT 



Boiml, 0, m? ) + 2m? i?o(m? , A^ m? 



Bo{m1^,ml, vr?) + (m? - - m?) BQim^.ml, vr?) 



q-ni gi q 



The four-squark interaction gives 



sS?(m^^ 



4a, 
3 vr 



' ^ - S®(mJ-J = ^sin4^, [Ao(mJ-J - A,{ml, 



m-gmq cos 26'g So(m?^, m?, m;^, 



where y4o(m^) is the standard one-point function in the convention of [|14 



(27) 



(2J 



(29) 
(30) 

(31) 



^The gluon loop due to the qqgg interaction gives no contribution because it is proportional to A^lnA 0. 
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3.3 Shift from bare to on— shell parameters 



We next fix the shift from the bare to the on-shell couplings ^O'f^^^ in Eq. (|TBp. We follow the 



procedure which was first given in |jTT|. From Eqs. (|T2]) and (P) we get for the squark decays 
into or {i = 1, 2) 



21 



COS 26 a 



6G! 



21 



with 5G/ obtained by varying Eqs. (^2|-|5^), e. g 

f Amt Sq 5mt 



S{mtAt) s„ - /iCa 6mt 
2m,„s^ I 5{mtAt)sa — iJiCaSmt irritSaSmt 



( 



6{mf,Ab)ca — fiSaSnib 



6G^ = : 



(32) 



(33) 
(34) 



with S/3 = sin/?, = cos/5, Sq, = sin a, Cq, = cos a, and a the Higgs mixing angle. 6Gi is 
obtained from Eqs. ( |33D and (0) by: 



SGf = (5Gl with a^a + Ti/2) . 



(35) 



For the couplings to the A^ boson we have explicitly 



5G 



He) 



213 2mw 

where cot P (tan/3) is for g = t (6). 

For the decay tj bj (/c = 4) we get 



[ 6{mgAq) {cot j3, tan /?} + Sniq ] , 



SG 



i{c) 
ij4 



7^* 6G4 in' 



b\T 



-iyGl>^,59i-{-iyGl^,,S9i, J^^, 



with 



6Ga 



V2 



g ( 2mbSmb tan P + 2mtSmt cot P 5 (TUbAb) tan P + fiSTUb 
I 6 {rrit At) cot P + ^5mt 2{6mtmb + mtdmb) / sm2P 



(36) 



(37) 



(38) 



and analogously the expression for bi tj H according to (|T 

6mg is the shift from the bare to the pole mass of the quark q and has two contributions 
(Fig. Id). The gluon exchange contribution is 



5o(m2,0,m2)-5i(m^,0 



r/2 



(39) 



and the gluino contribution is 



6mf 



[nig [Bi{ml,ml,mlJ + Bi{ml,ml,m^-^) 
+ m3sin26'q [5o(mJ,m?,m?J - 5o(mJ,m?,m|)] } . 



(40) 



The parameter r in Eq. (|39|) exhibits the dependence on the regular ization scheme. As r does 
not cancel in the final result, we have to use dimensional reduction |T2| (r = 0) which preserves 
supersymmetry at least at two-loop order. 



For the renormalization of rriqAq we define the on-shell parameters Mq fj ^ and ^ in terms 
of squark pole masses m^. and on-shell mixing angels 9q (which will be defined below) using 
the tree-level relations (|^) to (|^). We thus get |^ 



5{mgAq) = \ (5m? - 5m?) ^wildq + (m? 
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m^ 
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cos 26*^ 5^^^ + yU {cot /5, tan/5} 5mg. (41) 



where again cot /3 (tan/?) is for g = t (6). bmh, is given by 



Re 



sSfK) + sifK) + sirK) 



<9) 



(42) 



with 



^(9) 
Hi 



II y'^qi 



qi' 



2a 



m? {25o(m? , 0, m? ) + 5i(m? , 0, m? )}, 



3 TT 

'a TT 



{Ao(m^)+m?B,(m?,mi,m: 



m^mg 



sin2^g-] fio("^|,"^|,"^5)}, 



{cos'2% Ao(m? ) + sin2 2^,- Ao(m?,)}, i ^ i' . 



(43) 

(44) 
(45) 



We also have to renormalize the squark mixing angle. This problem was first solved in [|T^: 
there the counterterm 59q was fixed in the process e^e~ — q^q2 such that it cancels the off- 
diagonal part of the squark wave-function corrections. 66q = dOf'' + SOf^ is thus given by 



a.. 



sin 4:9a 



Git ml 



Ql 



ml 



q\' 



(46) 



i'g) 



m-gmq 



a. 



3vr lUml 



92 ' 



5o(m-2, mg, m^)vu- 5o(m-, , mg, mj V22 



91' 



(47) 



where Vij comes from the Z^q^q* couplings, Vn = 4(/|^ cos^ 9g — s^e^) and -522 = sin^ 9g — 
s^Cq) with = sin^^jy. We will use this scheme in what follows. There are also other 
possibilities of defining the on-shell squark mixing angle. In [IC], for instance, 69q was fixed 
such that the renormalized self-energy of the squarks remains diagonal on the qi mass shell. 
Similar conditions were used in [|^, 0. For comparison we list the counterterms 69q in these 
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schemes in terms of the squark self-energy S12: 



1 Re {wiiSi2(m|) - V22^i2{m^ 



66g([15]) = z 5 ^) (used in this paper) 



Re{Si2K)} 

ml - ml 

Re{Si2(g^)} 
ml - ml 



6e,i[17]) = I ^ '^^K (48) 

2 ml - ml 

The differences between them are ultraviolet finite. In Fig. |^ we compare 56^ of |TI], W, ^ 
to that of our scheme (||15[)- As can be seen, the numerical differences between the various 
schemes are very small (< 1%). 

There is yet another subtlety that has to be taken into account: At tree-level and in the DR 
renormalization scheme S\J{2)l symmetry requires that the parameter Mq in the t and b mass 
matrices have the same value. This is, however, not the case at loop-level in the on-shell 
scheme due to different shifts in the t and in the b sectors [|11|, [1^. In this paper we take 

— "^"^s^ 9g + ml sin^ 9g — m| cos 2/? {I^j^ — sin^ 9w) — as the on-shell parameter 



in the q sector (g = t,b). This leads to a shift of M? [[TT 



M? (6) = M?(t) + 5M? (t) - 5M?(6) (49) 

with 

SMq^q) = 5ml cos^ 9q + Sml sin^ 9g — (ml — ml) sin 29g S9g — 2mqSmq. (50) 

The underlying SI]{2)l symmetry is reflected in the fact that the shift 5M|(t) — 5M?(6) is 
finite. 



3.4 Real gluon emission 



In order to cancel the infrared divergence we include the emission of real (hard and soft) gluons 
(Fig. le): 



Jo + /i + m- Joo + m] III + X Jqi 



(51) 



Again, X = mf + — m\j . The phase space integrals /„, and J„m have {mi^mj^mn^) as 



arguments and are given in WA 



We have checked explicitly that the corrected decay width F of Eq. (|T6D is ultraviolet and 
infrared finite. 
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4 Numerical results and discussion 



Let us now turn to the numerical analysis. As input parameters we use m^_^ , m^^ , cos 9^, tan f3, 
H, rriA, and rrig. From these we calculate the values of the soft SUSY-breaking parameters 
MQ^i), MQib), Mfj f), and A,b according to Eqs. (|) to (|) together with Eqs. (|9D and (|0D, 



taking = 1.12 MQ(t) and Af, = At for simplicity to fix the sbottom sector. Moreover, we 
take mt = 175 GeV, = 5 GeV, = 91.2 GeV, sin^ 9w = 0.23, a(m^) = 1/129, and 
as(m^) = 0.12. For the running of we use as{Q) = 127r/[(33 — 2nj) ln((5^/A^^.)] with Uf 
the number of quark flavors. We take = as{mga) for the decay except for the shift in 
Eqs. (^) and (^) for which we take = as{MQ{t)). For the radiative corrections to the 
and masses and their mixing angle a {—^ < a < ^ hj convention) we use the formulae of 
IIT9IPI; for those to we follow |20]. In order to respect the experimental mass bounds from 



LEP2 |2T[ and Tevatron |2|] we impose m^o > 90 GeV, mf^-f,^ > 85 GeV, and rrig > 300 GeV. 
Moreover, we require 6p {t — b) < 0.0012 |^ from electroweak precision measurements using 
the one-loop formulae of H] and A^ < 3 (M|(t) + M? + m|^J, Al < 3 {M^{b) + M? + m^J 

with rrijj^ = {m\ + m|) cos^ ~\ ^1 aiid '^Hi = + ^1) sin^ 1^ ~\ [Hi fi'om tree-level 
vacuum stability. 

As a reference point we take m^-^ = 250 GeV, m^^ = 600 GeV, cos6'f = 0.26 (Qi ~ 75°), 
tan/3 = 3, = 550 GeV, = 150 GeV, and rrig = 600 GeV. This leads to mi = 564 
GeV g, m^^ = 627 GeV, cos^^ = 0.99, At^ = -243 GeV, niho = 100 GeV, m^o = 162 GeV, 
sin a = —0.58, and mH+ = 164 GeV. Thus t2 can decay into ti + {h^,H'^,A^), and 61,2 can 
decay into iiH~ . 

We first discuss the parameter dependence of the widths of ^2 decays into neutral Higgs bosons 
by varying one of the input parameters of the reference point. We define the SUSY-QCD 
corrections as the difference between the corrected width F of Eq. (^) and the tree-level width 
ro of Eq. (IID. 

Figure |^ shows the tree-level and the SUSY-QCD corrected widths of the decays ^2 — ^ + 
H'^, y4°) as a function of m^^. The relative corrections ST/T^ = (F — r°)/r° are about —10% 
for the decay into and —9% to —62% for the decay into A^. The corrections for t2 — > tiH^ 
are —9%, —45%, and +45% for m^^ = 420, 670, and 900 GeV, respectively. The spikes in the 
corrected decay widths for m^^ = 775 GeV are due to the ^2 — ^ tg threshold. The different 
shapes of the decay widths can be understood by the wide range of the parameters entering 
the Higgs couplings to stops. In the range m^^ = 300 GeV to 900 GeV, we have At = 144 GeV 
to -889 GeV and sin a = -0.52 to -0.73 {rriho = 81 GeV to 114 GeV, and mj^o = 163 GeV to 
170 GeV). 

Figure @ shows the cos 9t dependence of the tree-level and the corrected widths of ^2—^*^1 + 
y4°) decays. Again the shapes of the decay widths reflect their dependence on the 

^Notice that |20j have the opposite sign convention for the parameter /i. 

•^Notice that at tree-level one has m^^ — 560 GeV because Mq = 558 GeV for both the i and b mass matrices; 
at 0{as), however, one gets M^ii) = 558 GeV and Mq(6) = 563 GeV. 
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underlying SUSY parameters in a characteristic way. In particular we have At = 1033, 183, 
-666 GeV and sin a = -0.748, -0.565, -0.726 for cos^^- = -0.7, 0, 0.7, respectively. Apart 
from the points where the tree-level decay amplitudes vanish the relative corrections range 
from -40% to 20%. 

In Fig. ^ we show the tree-level and the SUSY-QCD corrected decay widths as a function of 
iriA. For ruA = 100, 200, 300 GeV we have m^o = 85, 104, 105 GeV, mno = 128, 207, 304 GeV, 
and sin a = —0.87, —0.45, —0.37, respectively. The corrections to T^(t2 — > tih^) range from 
-15% to -7% for TUA = 100 GeV to 400 GeV. Those to r°(t2 ^ iiH°) are -50% to -22% for 
niA > 114 GeV, and those to r°(t2 ^ hA^) are about -25%. 

As for the dependence on the gluino mass, as can be seen in Fig. ^, the gluino decouples very 
slowly: in the range m-g = 300 GeV to 1500 GeV 6T/T° varies from (-9%, -37%, -28%) to 
(-7%, -16%, -14%) for the decays ^ 1 1 + A°), apart from the ^2 tg threshold 

at rUg = 425 GeV. 

As for the dependence on tan/3, we get 1(^2 ^ ti/i°) = 2.68, 2.09, 1.42 GeV with ST/T^ ~ 
-10%, -7%, -5% for tan^g = 3, 10, 30, respectively. Likewise, we get F = 0.67, 1.61, 2.45 GeV 
with (5F/F0 ^ -27%, -19%, -17% for the decay into and F = 2.1, 2.64, 2.92 GeV with 
(5F/F° ^ -22%, -19%, -18% for the decay into A°, respectively. 

Let us now turn to the sbottom decays. We start again from the reference point given above. 
For the decay h iiH~ we get F = 3.88 GeV with 5F/F° = -24%, and for the decay 
h iiH- we get F = 0.08 GeV with 5F/F° = +87%. As in our examples the width of 
the latter decay is usually quite small (because 62 — &k and ti ~ tj^) we will discuss only the 
parameter dependence of the hi decay. 

Figure |^ shows the tree-level and the SUSY-QCD corrected widths of this decay as a function 
of mi^. The SUSY-QCD corrections are about —25%. Notice that at tree-level we have m^^ = 
556 GeV to 566 GeV for mi^ = 85 GeV to 400 GeV, whereas at C(a,) we have m^^ = 561 GeV 
to 570 GeV. Therefore, the thresholds at tree-level and one-loop level are slightly different. 
The dependence on the stop mixing angle is shown in Fig. |^ for tan/3 = 3 and 10, and the 
other parameters as given above. (For jcos^l > 0.72 the decay bi —>■ tiH^ is kinematically 
not allowed.) In case of tan/5 = 3, the SUSY-QCD corrections range from about —40% to 
26%, with 5T /T^ > for cos^^f < — 0.6. In case of tan/3 = 10 ^F/F" is much larger. For 
cosOf > 0.5 and tan/? = 10 we even get a negative corrected decay width. This is mainly due 
to a large contribution stemming from the term S{mbAb) ~ /itan/35mb of Eq. (0) and was 
already mentioned in [|ll|]. The same problem can occur for the decays 62 &i + {h^, H^, A^) 



and ^2 biH'^ which may be important for large tan/3 (due to the large bottom Yukawa 
coupling and the large 61-62 mass splitting). 

We have also studied the dependence on niA- In the case tan/? = 3 (10) we have found that 

5F/F°(6i ^ iiH-) 20% (-40%) for niA = 100 GeV to 285 GeV and the other parameters 

as given above. 

As for the dependence on the gluino mass, 5F/F'^(6i —>■ tiH~) ranges from —26% to —14% 
(-47% to -39%) for m-g = 300 GeV to 1500 GeV and tan/3 = 3 (10). 

In [|r^ a numerical analysis for the decays ^2 ~* ^1 + {h^, A^) was made. Whereas we agree 
with their figure for t2 —>■ tiA^, we find a difference of about 10% in both the tree-level and 
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the corrected widths of t2 —>■ tih^. This may be due to a different treatment of the radiative 
corrections to the mass and mixing angle a. 



5 Conclusions 



We have calculated the 0{as) SUSY-QCD corrections to the decay widths of q2 ^ Qi + 
y4°) and ^j-^^ (•? ~ ^) on-shell scheme. We have taken into account 

appropriate shifts for the soft SUSY-breaking parameters defined in terms of on-shell squark 
masses and mixing angles. It has turned out that the SUSY-QCD corrections are mostly 
negative and of the order of a few ten percent and should therefore be taken into account. 
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Appendix: Squark — Squark — Higgs Couplings 



squark - squark - 



— ^ T ^, 



\Pi TTlq hq 



So 
-Co 



gm 



(52) 



.a, 



for I Jq^j^ } type squarks respectively. We use the abbreviations c^y = cos6w, Sq, = sin( 
Cq, = cos a, Sa+/3 = sin(a+/3), Cj^ = Ili—eq sin^ 6^r, and = Cq sin^ O^r. a is the mixing angle 
in the CP even neutral Higgs boson sector, hq are the Yukawa couphngs: 



h. 



grrib 



my/ sin /? 



mw cos/3 



(53) 



*We thank A. Djouadi and W. HoUik for correspondence on this point. 
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(b) squark - squark - 



(54) 



(c) squark - squark - 



2m\Y 



cot /3 
tan/3 



_^ /cot/3 



« \taii/3 




'1 



(55) 



(d) squark - squark - H"^ 



V2r 



ml tan P + mf cot /3 — sin 2/3 m;, (>!;, tan (3 -\- ji) 

mt {At cot /3 + A*) 



2mf TTtj, 
sin 2/3 



(56) 
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/ ^3 



(a) 



/ ^3 




0° V M 



(b) 




Figure 1: Fcynman diagrams relevant for the 0{as) SUSY-QCD corrections to squark decays 
into Higgs bosons. 
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Figure 2: Differences in ScosOf between [|I 



for m 



m 



250 GeV, 



16, 17[ and our sclieme |1T5| 
600 GeV, and m-g = 600 GeV. For 5cos^t(Q2)([i6]) we liave taken 



t2' 
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mi, [GeV] 

Figure 3: Tree-level (dashed lines) and 0(as) SUSY-QCD corrected (solid lines) decay widths 
of i2 h + {h°, H^, A^) as a function of m^-^, for rrii^ = 250 GeV, cos% = 0.26, fi = 550 GeV, 
tan/3 = 3, ttia = 150 GeV, and rrig = 600 GeV. The grey area is excluded by the bound 
ruho > 90 GeV. 
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Figure 4: Tree-level (dashed lines) and 0{as) SUSY-QCD corrected (solid lines) decay widths 
of t2 ^ ti + (/i° H^, as a function of cos ^f, for m^^ = 250 GeV, m^^ = 600 GeV, = 550 GeV, 
tan/3 = 3, rriA = 150 GeV, and rrig = 600 GeV. The grey areas are excluded by the constraints 
given in the text: 5p{t-h) > 0.0012 in (a), rriho < 90 GeV in (b), and unstable vacuum in (c). 




Figure 5: Tree-level (dashed lines) and 0{as) SUSY-QCD corrected (solid lines) decay widths 
of t2 ii + {h^, H^, A^) as a function of tua, for m^-^ = 250 GeV, m^-^ = 600 GeV, cos^^- = 0.26, 
II = 550 GeV, tan/3 = 3, and nig = 600 GeV. The grey area is excluded by the bound 
rriho > 90 GeV. 



16 










1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

tihy 




-10 
















-20 












o 

l-H 


-30 




l-H 










-40 








-50 


-i 1 1 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 ~ 



400 600 800 1000 1200 1400 

nig [GeV] 



Figure 6: SUSY-QCD corrections (in %) to the widths of ^ h + {h^, H^, as a 

function of rrig, for m^^ = 250 GeV, m^^ = 600 GeV, cos6'j = 0.26, /i = 550 GeV, tan/3 = 3, and 
ruA = 150 GeV. 




100 200 300 400 

mi, [GeV] 

Figure 7: Tree-level (dashed line) and 0{as) SUSY-QCD corrected (solid line) decay widths 
of bi tiH~ as a function of m^^, for m^-^ = 600 GeV, cos6'j = 0.26, /x = 550 GeV, tan/3 = 3, 
tua = 150 GeV, and rrig = 600 GeV. The insert zooms on the different thresholds at tree- and 
one-loop level. 
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Figure 8: Tree-level (dashed lines) and 0{as) SUSY-QCD corrected (solid lines) decay widths 
of bi tiH^ as a function of cosO^, for m^^ = 250 GeV, m^^ = 600 GeV, /i = 550 GeV, 
VIA = 150 GeV, rrig = 600 GeV, and tan/3 = 3, 10. The grey area is excluded for tan (3 = 3 by 
the bound m^o > 90 GeV. 
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